Calculation of binary magnetic properties and potential energy curve in xenon dimer: second virial coefficient of (129)Xe nuclear shielding

Quantum chemical calculations of the nuclear shielding tensor, the nuclear quadrupole coupling tensor, and the spin-rotation tensor are reported for the Xe dimer using ab initio quantum chemical methods. The binary chemical shift delta, the anisotropy of the shielding tensor Delta sigma, the nuclear quadrupole coupling tensor component along the internuclear axis chi( parallel ), and the spin-rotation constant C( perpendicular ) are presented as a function of internuclear distance. The basis set superposition error is approximately corrected for by using the counterpoise correction (CP) method. Electron correlation effects are systematically studied via the Hartree-Fock, complete active space self-consistent field, second-order M?ller-Plesset many-body perturbation, and coupled-cluster singles and doubles (CCSD) theories, the last one without and with noniterative triples, at the nonrelativistic all-electron level. We also report a high-quality theoretical interatomic potential for the Xe dimer, gained using the relativistic effective potential/core polarization potential scheme. These calculations used valence basis set of cc-pVQZ quality supplemented with a set of midbond functions. The second virial coefficient of Xe nuclear shielding, which is probably the experimentally best-characterized intermolecular interaction effect in nuclear magnetic resonance spectroscopy, is computed as a function of temperature, and compared to experiment and earlier theoretical results. The best results for the second virial coefficient, obtained using the CCSD(CP) binary chemical shift curve and either our best theoretical potential or the empirical potentials from the literature, are in good agreement with experiment. Zero-point vibrational corrections of delta, Delta sigma, chi (parallel), and C (perpendicular) in the nu=0, J=0 rovibrational ground state of the xenon dimer are also reported.     

MM-PBSA free energy analysis of endo-1,4-xylanase II (XynII)-substrate complexes: binding of the reactive sugar in a skew boat and chair conformation

The binding of xylotetraose in different conformations to the active site of endo-1,4-beta-xylanase II (XynII) from Trichoderma reesei was studied using molecular dynamics (MD) simulations and free energy analyses employing the MM-PBSA (Molecular Mechanics-Poisson-Boltzmann Surface Area) method. MD simulations of 1 ns were done for the substrate xylotetraose having the reactive sugar, which is bound in the -1 subsite of XynII in the 4C1 (chair) and 2So (skew boat) ground state conformations, and for the transition state of the XynII catalysed hydrolysis of the beta-glycosidic linkage. According to the simulations and free energy analysis, XynII binds the substrate with the -1 sugar in the 2So conformation 59.8 kJ mol(-1) tighter than the substrate with the sugar in the 4C1 conformation. The reactive 2So conformation resembles closely the reaction transition state and has the breaking glycosidic bond in a pseudo-axial orientation ready for facile bond cleavage. The transition state was calculated to be bound 77.1 kJ mol(-1) tighter than the 4C1 ground state conformation. The molecular mechanical interaction energy between the enzyme and the reactive pyranoside unit at the -1 subsite was 75.7 kJ mol(-1) more favorable for the binding of the 2So conformation than the 2C1 conformation, explaining the clearly tighter binding of the reactive structure The results of this study indicate that in the Michaelis complex XynII, a member of the family 11 xylanase, the substrate is bound in a skew boat conformation and in the catalytic reaction, the -1 sugar proceeds from the 4C1 conformation through 2So to the transition state with the -1 sugar in the 2,5B conformation.     

Cytochrome c–dimyristoylphosphatidylglycerol interactions studied by asymmetrical flow field-flow fractionation

Lipid membranes are well recognized ligands that bind peripheral and integral proteins in a specific manner and regulate their function. Cytochrome c (cyt c) is one of the partner peripheral protein that binds to the lipid membranes via electrostatic and hydrophobic interactions. In this study, asymmetrical flow field-flow fractionation (AsFlFFF) was used to compare the interactions of cyt c with the acidic phospholipid 1,2-dimyristoyl-sn-glycero-3-phospho-rac-glycerol (DMPG), oleic acid (OA), and sodium dodecyl sulfate (SDS). The influence of pH and the cyt c–lipid molar mass ratios were evaluated by monitoring the diffusion coefficients and particle diameter distributions obtained for the free and lipid-bound protein. The hydrodynamic particle diameter of cyt c (pI 10) was 4.1 nm at pH 11.4 and around 4.2 nm at pH 7.0 and 8.0. Standard molar mass marker proteins were used for calibration to obtain the molar masses of free cyt c and its complexes with lipids. AsFlFFF revealed the binding of cyt c to DMPG and to OA to be mainly electrostatic. In the absence of electrostatic interactions, minor complex formation occurred, possibly due to the extended lipid anchorage involving the hydrophobic cavity of cyt c and the hydrocarbon chains of DMPG or SDS. The possibility of the formation of the molten globule state of cyt c, induced by the interaction between cyt c and lipids, is discussed.Electronic Supplementary Material Supplementary material is available for this article at     

Semiempirical PM5 molecular orbital study on chlorophylls and bacteriochlorophylls: comparison of semiempirical, ab initio, and density functional results

The semiempirical PM5 method has been used to calculate fully optimized structures of magnesium-bacteriochlorin, magnesium-chlorin, magnesium-porphin, mesochlorophyll a, chlorophylls a, b, c(1), c(2), c(3), and d, and bacteriochlorophylls a, b, c, d, e, f, g, and h with all homologous structures. Hartree-Fock/6-31G* ab initio and density functional B3LYP/6-31G* methods were used to optimize structures of methyl chlorophyllide a, chlorophyll c(1), and methyl bacteriochlorophyllides a and c for comparison. Spectroscopic transition energies of the chromophores and their 1:1 or 1:2 solvent complexes were calculated with the Zindo/S CIS method. The self-consistent reaction field model was used to estimate solvent shifts. The PM5 calculations predict planar structure of the porphyrin ring and central position of the four coordinated magnesium atoms in all pigments studied, in accord with the experimental, ab initio, and density functional results, a significant improvement as compared to the older semiempirical PM3 approach. Only small differences in PM5 and B3LYP/6-31G* or Hartree-Fock/6-31G* minimum energy geometries of the reference molecules were observed. Calculations show that in 1:1 solvent complexes, where the magnesium atom is five coordinated, the magnesium atom is shifted out of the plane of the porphyrin ring towards the solvent molecule, while the hexa coordinated 1:2 complexes are again planar. The PM5 method gives atomic charges that are comparable with those obtained from the Hartree-Fock/6-31G* and B3LYP/6-31G* calculations. The single point ZINDO/S CIS calculations with PM5 minimum energy structure gave excellent correlations between calculated and experimental transition energies of the chlorophylls and bacteriochlorophylls studied. Such correlations may be used for prediction of transition energies of the chromophores in protein binding sites. Calculations also predict existence of dark electronic states below the main Soret absorption band in all chromophores studied. The results suggest that the semiempirical PM5 method is a fairly reliable and computationally efficient method in predicting molecular parameters of porphyrin-like molecules.     

Inclusive ϱ0,K *±(892) andf production in\bar p p Interactions at 12 GeV/c

Inclusive production of ?⁰,K ^*±(892), andf is studied in \(\bar p\) p interactions at 12 GeV/c. The inclusive cross sections for ?⁰,K ^*±(892), andf are found to be 6.7±0.3 mb, 1.0±0.2 mb, and 1.4±0.3 mb, respectively. The differential cross sections are presented as a function of c.m. rapidity, Feynmanx and square of the transverse momentump _T ² . Comparison with the correspondingpp data shows some interesting differences which can be attributed to the \(\bar p\) p annihilation. The results are compared with the predictions of the quark fusion model.     

Supramolecular chirogenesis in zinc porphyrins: interaction with bidentate ligands,formation of tweezer structures,and the origin of enhanced optical activity

The complexation behavior, binding properties, and spectral parameters of supramolecular chirality induction in the achiral host molecule, syn (face-to-face conformation) ethane-bridged bis(zinc porphyrin), upon interaction with chiral bidentate guests (diamines and amino alcohols) have been studied by means of UV-vis, CD, fluorescence, (1)H NMR, and ESI MS techniques. It was found that the guest structure plays a decisive role in the chirogenesis pathway. The majority of bidentate ligands (except those geometrically unsuitable) exhibit two major equilibria steps: the first guest ligation leading to formation of the 1:1 host-guest tweezer structure (K(1)) and the second guest molecule ligation (K(2)) forming the anti bis-ligated species (1:2). The second ligation is much weaker (K(1) > K(2)) due to the optimal geometry and stability of the 1:1 tweezer complex. The enhanced conformational stability of the tweezer complex ensures an efficient chirality transfer from the chiral guest to the achiral host, consequently inducing a remarkably high optical activity in the bis-porphyrin.     

On applicability of PCA,voxel-wise variance normalization and dimensionality assumptions for sliding temporal window sICA in resting-state fMRI

Subject-level resting-state fMRI (RS-fMRI) spatial independent component analysis (sICA) may provide new ways to analyze the data when performed in the sliding time window. However, whether principal component analysis (PCA) and voxel-wise variance normalization (VN) are applicable pre-processing procedures in the sliding-window context, as they are for regular sICA, has not been addressed so far. Also model order selection requires further studies concerning sliding-window sICA. In this paper we have addressed these concerns. First, we compared PCA-retained subspaces concerning overlapping parts of consecutive temporal windows to answer whether in-window PCA and VN can confound comparisons between sICA analyses in consecutive windows. Second, we compared the PCA subspaces between windowed and full data to assess expected comparability between windowed and full-data sICA results. Third, temporal evolution of dimensionality estimates in RS-fMRI data sets was monitored to identify potential challenges in model order selection in a sliding-window sICA context. Our results illustrate that in-window VN can be safely used, in-window PCA is applicable with most window widths and that comparisons between windowed and full data should not be performed from a subspace similarity point of view. In addition, our studies on dimensionality estimates demonstrated that there are sustained, periodic and very case-specific changes in signal-to-noise ratio within RS-fMRI data sets. Consequently, dimensionality estimation is needed for well-founded model order determination in the sliding-window case. The observed periodic changes correspond to a frequency band of ≤ 0.1 Hz, which is commonly associated with brain activity in RS-fMRI and become on average most pronounced at window widths of 80 and 60 time points (144 and 108 s, respectively). Wider windows provided only slightly better comparability between consecutive windows, and 60 time point or shorter windows also provided the best comparability with full-data results. Further studies are needed to determine the cause for dimensionality variations.